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An algorithm for the numerical solution of the Schrodinger equation in the case of a time de¬ 
pendent potential is proposed. Our simple modification upgrades the well known method of Koonin 
while negligibly increasing the computing time. In the presented test the accuracy is enhanced by 
up to an order of magnitude. 
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The microscopic description of many-body systems like atoms or nuclei is based on a many-body Hamiltonian. 
The related wave functions are given by Slater determinants for fermions. In the case of time dependent processes like 
atomic or nuclear collisions, nuclear fission or fusion, however, the situation is too complex due to the great number of 
degrees of freedom. Therefore, in most cases, a collective coordinate is introduced according to the essential physical 
properties of the considered system. This procedure leads to a macroscopic model with one degree of freedom, which 
is governed by an effective one-body Schrodinger equation including a time dependent potential in the considered 
examples [Q . The well known coordinate representation reads 
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A momentum dependent potential V, leading to an integro-differential equation, is excluded here. In almost all 
practical cases eq. (]^) has to be evaluated numerically. The standard technique essentially consists of the following 
two steps: first apply a suitable scheme for the space discretization and then perform the time integration. 

The algorithm presented here is only related to the second step: The time integration in the case of an explicitly 
time dependent potential H(r,t). Our proposed convenient modification upgrades the standard method to a much 
more efficient version while negligibly increasing the computing time. 

Our experience is due to the description of ternary fission, i.e. a fission process accompanied by the emission of an 
a-particle, with results to be published elsewhere. However, the presented algorithm might be of broader interest, 
not necessarily restricted to this nuclear physics theme. 

In this letter we will discuss the application for the case of a cylindrically symmetric, time dependent potential, 
essentially following Koonin’s procedure [P, and present a test using an analytically solvable example. 

As usual, in a cylindrically symmetric case, the angle dependence on ip of the wave function is separated as 
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and a grid is defined by 


= j = l,...,Np, 

Zk = kAz, k =-N^,... ,Nz . 
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For simplification we introduce now 


9j,k — y/Pj '4’j,k 
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leading to the space-discretized equation 

d 

^ ^ ^^S)],k = {vg)j^k + {hg)j,k , (5) 

where the Hamiltonian H is split in a “vertical” {v) and a “horizontal” (h) part 
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with the abbreviations 


{vg)j,k = - 


{hg)j,k = - 


^3 9j+3,k 2 gj^k “t“ Cj —1 gj — i^k 
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Ap' 


Uj^k gj,k 


gj,k+i — 2 gj,fc + gj,k-i 
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Az 2 


Uj^k gj,k 
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A,k - ^Vj,k + 
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Up to this point no changes are made in comparison to Koonin’s algorithm. We solve eq. (^) iteratively, with times 
tn = nAt, using a Taylor expansion of gjkitn) = gj^\ up to and including (At)2 and obtain 


= (^ 1 - 1 i? (At)i 

- ^ _^y ^ 3^^ gin) ^ (8) 


where all terms on the right hand side are evaluated at time tn- The time derivative Vj^k = dVj^k/dt is included to 
ensure consistency up to second order contributions. 

Inserting on the right hand side is not an advisable method to determine because of serious numerical 

instabilities. Therefore, an “alternating direction implicit method” j^] is used, leading to a modified version of 
Koonin’s algorithm 

(l + ^AAi + ^F(A()>) 

where again all terms on the right hand side are evaluated at time tn- Our result may be verified by an expansion in 
terms of At up to second order. There are no restrictions on the commutation of [u, h]. By neglecting terms dependent 
on Vj^k, our expression reduces to the one used by Koonin et al. in [^. 

For numerical convenience we now define by 


leading to 
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In order to determine , and thereafter , one has to invert tridiagonal matrices only, which can be performed 

by a suitable algorithm Q| . 

A central point for the efficiency of our modification is that the determination of 
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up to order (At)° poses no appreciable effort. Anyway, the values of Vj^k have to be calculated for all time steps in 
order to determine v and h. It is deemed unsuitable to compensate this by determining v and h at an intermediate 
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time step t + At/2 because one has to calculate the potential at this additional time step, as well. In many cases this 
produces appreciable effort. For example, in our calculations concerning the ternary fission it took approximately the 
same time to calculate for one value of n on the whole grid as to perform one iteration step from to . 

To test our algorithm, with regard to the modification due to keeping terms dependent on we consider the 
following, explicitly time dependent potential 

V (r, t) = i r ^ - 2nu;^ t . (13) 

At time t — 0 this is the potential of the harmonic oszillator with frequency u. For the initial wave function we choose 
an eigenfunction of this harmonic oszillator 

</(r,t = 0) = exp ^ , (14) 

where /? = Mu/h, i S {0,1,2,...}, Hi = Hermite polynomial, and ci = constant of normalization. 

The analytical solution is given by 

(r, t) = (/(r, t = 0) exp (—[£ + 3/2] t + . (15) 

The results are displayed in Table || for different radial quantum numbers i (different numbers of nodes of the wave 
function). Both parts of the wave function—real and imaginary—are considered separately. We use a; = 2, At = 1/60 
with 60 time steps and grid parameters Np = = 64 and Ap = Az = 0.25. The enhancement in accuracy becomes 

as large as an order of magnitude. The norm of the wave function is (nearly) conserved, since the algorithm—modified 
or not—is nearly unitary |Q. 

We note that the considerations for a one dimensional time dependent potential are quite the same. One obtains 

(a;, t + At) = ^ 1 + ^ Ft At + ^ 1/ {Atf ^ 

(^l-±HAt-^V{Atf^cl,ix,t) + 0 (At3) . 

(16) 

Again the calculation of V {x,t) poses nearly no effort. In contrast to the ordinary method, where the V terms are 
neglected |Q|, the above expression is accurate up to order (At)^. The approximation ( ^ is unitary; it automatically 
maintains the normalization of </>. We note that no alternating direction method is needed in the one dimensional 
case. 

We summarize: In the case of a time dependent potential the standard algorithm of Koonin can easily be modified 
to obtain a greater accuracy. Therefore, terms depending on V have to be included to get a consistent time expansion 
up to order (At)^. The related algorithm has been shown to be superior in the presented numerical test. 
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TABLE I. Results of the numerical test. The relative error of the real and imaginary part of the wave function is displayed 
for different numbers of nodes. The modihcation of the Koonin’s standard algorithm lead to an enhancement in the accuracy 
of up to an order of magnitude. 


nodes 

standard 

Relative error of 
the real part: 

modihed 

standard 

Relative error of 
the imaginary part: 

modihed 

£ = 0 

10% 


0.4% 

4% 


0.2% 


10% 


0.6% 

4% 


0.4% 

i = 2 

1% 


0.3% 

40% 


5% 

£ = 3 

20% 


4% 

1.5% 


0.6% 
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